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Abstract —In this paper we examine a number of deployment 
issues which arise from practical considerations in massive 
multiple-input-multiple-output (MIMO) systems. We show both 
spatial correlation and line-of-sight (LOS) introduce an interfer¬ 
ence component to the system which causes non-orthogonality 
between user channels. Distributing the antennas into multiple 
clusters is shown to reduce spatial correlation and improve 
performance. Furthermore, due to its ability to minimize inter¬ 
ference, zero forcing (ZF) precoding performs well in massive 
MIMO systems compared to matched filter (MF) precoding 
which suffers large penalties. However, the noise component in 
the ZF signal-to-noise-ratio (SNR) increases significantly in the 
case of imperfect transmit channel state information (CSI). 

I. Introduction 

Global mobile data traffic is expected to increase at a 
compound annual growth rate of 61% from 2013-2018 (T). 
To accommodate these demands, engineers are proposing new 
technologies for enhanced network capabilities. A symbio¬ 
sis of three key technologies: multiple-input-multiple-output 
(MIMO), millimeter wave (mmWave) and small cell technol¬ 
ogy, has been proposed 10 , 0 in the preliminary research to 
fifth generation wireless system (5G) standardization. 

Massive MIMO, which scales up the number of base station 
(BS) antennas by at least an order of magnitude, provides 
significant improvements in spectral efficiency, interference 
mitigation, data rate and robustness |4), 0 , relative to conven¬ 
tional MIMO. Additionally, as the number of antennas grows 
|| 6 l , user channels start to become mutually orthogonal. As a 
result, massive MIMO can be implemented with inexpensive, 
low-powered amplifiers 0 . However, the deployment of such 
a large number of antennas introduces many practical limiting 
factors. 

This paper is intended to highlight several deployment 
issues for downlink (DL) massive MIMO by examining prac¬ 
tical scenarios. We consider the effects of spatially correlated 
antenna topologies, imperfect channel state information (CSI), 
distributed antenna systems and line-of-sight (LOS) propa¬ 
gation on linear precoding performance in massive MIMO 
systems. All results are presented in terms of matched filter 
(MF) expected per-user signal-to-interference-plus-noise-ratio 
(SINR) and zero forcing (ZF) expected per-user signal-to- 
noise-ratio (SNR). 

The contributions of this paper are as follows: 

• We explore the impact of different antenna topologies on 
antenna array spatial correlation and system performance. 

• We demonstrate that distributing an antenna array into 
multiple clusters reduces spatial correlation and positively 


impacts system performance. 

• We consider a Rician fading channel to model the effects 
of LOS propagation, which we show to adversely affect 
massive MIMO system performance. 


II. System Model 
A. System Description 

We consider a single-cell massive multi-user (MU)-MIMO 
DL system with a total of M transmit antennas divided equally 
among N antenna clusters, jointly serving a total of K « M 
single-antenna users. At each antenna cluster, the antennas 
are assumed to be arranged as pairs of cross-polarized 
(x-pol) antennas. We assume time division duplex (TDD) 
operation with uplink (UL) pilots enabling the transmitter to 
estimate the DL channel via channel reciprocity. On the DL, 
the K single antenna terminals collectively receive the K x 1 
vector 

y = v / pG Tx + n , (1) 


where p is the transmit SNR, x is an M x 1 precoded data 
vector and n is a K x 1 noise vector with independent and 
identically distributed (i.i.d.) CAf(0 ,1) entries. The transmit 
power is normalized, E [||x|| 2 ] = 1, i.e., each antenna trans¬ 
mits at a power of -jA. The M x K UL channel matrix, G, is 
given by 


G = 


A’lRfHr,! 


Pn,1^ Hjv,! 


R?h 1>k 


Pn,k- R t 2 Hjv,k 


( 2 ) 


where H nt k £ C^ xl is the i.i.d. channel vector between the 
??th antenna cluster and the kth user, /?„ & is the link gain 
coefficient, modeling large-scale effects from antenna cluster 
n to user k, while R t is the spatial correlation matrix at each 
antenna cluster, assumed equal for all antenna clusters. 


B. Link Gain Model 

We consider the scenario where users are dropped at ran¬ 
dom locations within a circular coverage region. Antenna 
clusters are positioned equidistant on the periphery of the 
coverage region for N > 1 while for co-located antenna 
systems, all antenna elements are located in the centre of the 
circular coverage region. Path loss is then calculated from 
a conventional distance based model with i.i.d. log-normal 
shadowing: ALd ~ 7 , where L is i.i.d. log-normal shadowing 
with standard deviation (SD) er, d is the link distance, 7 is the 




path-loss exponent and A normalizes the 7,, ^ values so that 
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C. Spatial Correlation Model z - 

We consider a cross-polarized (x-pol) antenna configuration, 
with the spatial correlation matrix modeled via 0 

Ay 1 

l 



1 * 

R, = X pol 0 R, (3) dl I 
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where the MxM matrix R is the co-polarized (co-pol) spatial 
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correlation matrix, 0 represents the Hadamard product and 
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Xp 0 i is the MxM x-pol matrix given by 







Xp ol — 1 M_ 0 
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where 1 m is a 4r x 4/- matrix of ones, <5 denotes the cross- 

2 Z Z 

correlation between the two antenna elements in the x-pol 
configuration and (g> represents the Rronecker product. 


(a) (b) 

Fig. 1: (a) URA. (b) Cylindrical Array. 


normalizes the average power of the MF precoder. The ex¬ 
pected value of the ith user’s MF SINR is then given by lfl~2l 


D. Imperfect CSI Model 

We model CSI imperfections via an estimated channel 
matrix, G, given by (9) 

G = + y/l-^E, (5) 
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+ 1 


where E is independent and statistically identical to G and £, 
0 < £ < 1 . controls the accuracy of the CSI. 


E. Linear Precoders 


Due to their simplicity and optimality in massive MU- 
MIMO systems HOI, we examine the performance MF and 
ZF linear precoding techniques. 

1) MF Precoding: The MF precoder (also known as a 
matched beamformer (MBF) and maximum ratio transmission 
(MRT)) is the most computationally inexpensive precoding 
technique, allowing the design of many inexpensive antennas 
ideal for massive MIMO systems. MF precoding in a MU- 
MIMO system aims at maximizing the received power at 
each user while neglecting the effects of interference to the 
other co-scheduled users. As the number of transmit antennas 
increases without bound, for a fixed number of users and 
perfect CSI, MF precoding benefits from the law of large 
numbers, effectively eliminating all inter-user interference 
HD- As a consequence, for an infinite number of antennas, 
all user channels become mutually orthogonal and information 
capacity is maximized. 

For a MF precoder, the M x 1 precoded data vector, with 
CSI inaccuracies, is given by G2 


X = 



( 6 ) 


where q is the K x 1 data symbol vector, with E [||q|| 2 ] = 1, 
and 


7mf 


tr(G T G*) 

K 


(7) 


where P, = E [e*ej], e, is the ith column of E, g, is the ith 
column of G and the noise power is normalized to 1. 

2) ZF Precoding: The ZF precoding technique forces all 
intra-cell interference to zero by using coherent superposition 
of wavefronts to send null vectors to all other users. The 
ZF precoder is more computationally expensive than the MF 
precoder as it requires a matrix inverse of G T G*. However 
with large antenna numbers W = (G T G *)/M tends to the 
identity matrix, 1 ^, and the computation of the matrix inverse 
becomes trivial. 

For a ZF precoder, the M x 1 precoded data vector, with 
CSI inaccuracies, is given by 0 

x = - ? LG*(G T G*)- 1 q, (9) 


where 


7zf 
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K 


( 10 ) 


normalizes the average power of the ZF precoder. The ex¬ 
pected value of the ith users ZF SNR is then given by 


E [SNR,] 


-JI—p 2 
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where the noise power is normalized to 1 . 


III. Antenna Array Topologies 

It was shown in 0 that the additional benefits expected of 
massive MIMO systems are realized only when the number 
of antennas is on the order of a thousand. With such stringent 
requirements, antenna arrays can become unfeasibly large. 
























e.g., 256 co-polarized antenna elements positioned at a half¬ 
wavelength spacing at / = 2.6 GHz requires a uniform 
linear array (ULA) of 14.8 m. To combat space requirements, 
antenna arrays need to utilize both azimuth and zenith di¬ 
mensions m. We explore the impact on linear precoding 
performance of deploying antennas in uniform rectangular 
array (URA) and cylindrical antenna topologies, shown in 
Figures [Ta| and 1 1 hi respectively. 

We observe the effects of spatial correlation, we constrain 
the array dimension to l (the width and the height of URA 
topology). The cylindrical array can then be obtained by 
wrapping the URA around a virtual cylinder. Thus, in Figure 
US the circumference and radius of the uniform circular 
array (UCA) on the x,y -plane is equal to l and r = 2n/l, 
respectively. Spatial correlation for each antenna topology is 
then generated via 

R = AF 0 AF» (12) 


where AF rJ and AF@ are the azimuth and zenith domain 
antenna array factors of size P x 1 and Q x 1 respectively, 
given by Cl 


AFrf, = 


l e -i$i(A0,AS) 


AFa = 


1 e-ieHAG. 


e -j0Q-i(A0) 


(13) 

(14) 


AO) and 0 g (A 0) are the azimuth and zenith domain 
phase shifts of the pth and qth antenna’s angle of departure 
(AOD) respectively, with respect to a reference antenna, given 
by 


$p(A0, A 6) = kd p cos (<j) + A(f>) sin(0 + A 6), (15) 

O q (A0) = kd q cos(6 + Ad), (16) 

where k is the wavenumber, d p is the distance between the 
pth antenna and the reference antenna element in the azimuth 
domain, d q is the distance between the r/th antenna and the 
reference antenna element in the zenith domain, A</> is the 
azimuth domain AOD offset relative to its mean , and AO is 
the zenith domain AOD offset relative to its mean 0. We model 
the probability density functions (pdfs) of LOS and non-line- 
of-sight (NLOS) AODs as described by the 3rd Generation 
Partnership Project (3GPP) IIT51 3D channel model. 

All results are generated by simulating MF E[SINR,] and 
ZF E[SNRj], in © and ( II 11 . respectively, over many channel 
realizations. System parameters given in Table [I] where A is 
the wavelength. 

In Figures [2] and [3] we show the impact of URA and 
cylindrical antenna topologies on the MF and ZF precoding. 
We plot the cumulative density functions (CDFs) of per user 
SINR and per user SNR for MF and ZF, respectively. In 
each case of precoding technique, it can be observed that 
the URA topology has a marginally better performance than 
the cylindrical topology. This is a result of smaller inter¬ 
element distances in the azimuth domain of the cylindrical 
array topology compared with the URA. Note, that Maxwell’s 


Parameter 

Value 

Transmit SNR, p 

10 dB 

Log-normal shadowing SD, cr 

8 dB 

Path-loss exponent, 7 

4 

Link distance, d 

50-1000 m 

Maximum link gain coefficient, /3 max 

25 dB 

System frequency, / 

2.6 GHz 

X-pol parameter, 6 

0 .1 2 

Array dimension, l 

2A 

Azimuth domain AOD pdf, pa 0 (A</>) 

Wrapped Gaussian 

Zenith domain AOD pdf, pAe(AO) 

Laplacian 

NLOS mean azimuth AOD, (f) 

74.13° 

NLOS mean zenith AOD, 0 

18.20° 

NLOS azimuth AOD SD, a/^ 0 

1.29° 

NLOS zenith AOD SD, a a 9 

1.45° 

LOS mean azimuth AOD, 

64.57° 

LOS mean zenith AOD, 0 

8.91° 

LOS azimuth AOD SD, a a </> 

1.58° 

LOS zenith AOD SD, a ah 

1.45° 


TABLE I: System Parameters 



Fig. 2: MF E[STNR, ; ] CDF as a function of CSI accuracy, £, 
and antenna topology for M = 256 and K = 8 . 


equations fundamentally limit the effectiveness of 3D arrays 
such that only antennas on the surface of the array contribute 
to the information capacity 0. The large tails of the ZF 
precoder CDFs, in Figure [3| are due to the sub-optimality of 
the precoder at low SNR, where noise becomes dominant. The 
ZF precoder is more robust to the effects of spatial correlation 
than the MF precoder, i.e., the reduction in ZF SNR is much 
less than the reduction in MF SINR when spatial correlation is 
introduced. This is because spatial correlation acts as a form 
of interference to the system which the ZF precoder is able 
to minimize. The MF precoder, on the other hand, is more 
robust to the effects of imperfect CSI. The additional noise 
component in the denominator of the ZF SNR, in ( flTb . is 
increased significantly as £ is reduced. 

















and antenna topology for M = 256 and K = 8. 



antenna topology for M = 256 and K = 8. 

In Figures [4] and [5] we present the impact of varying angle 
spread (AS) on URA and cylindrical antenna topologies. We 
consider the extreme cases of 0.125x and 8x the 3GPP angle 
spreads (given in Tabled on MF and ZF precoding techniques. 
We observe negligible difference in ZF precoding performance 
for all angle spreads considered. However, the performance of 
MF precoding is improved slightly with a larger angle spread, 
e.g., the cylindrical antenna topology shows an increase of 
median value by approximately 2 dB E[STNR.;]. This is a result 
of a larger angle spread increasing the number of independent 
paths from each antenna, thus decreasing spatial correlation. 

IV. CO-LOCATED VS DISTRIBUTED ANTENNA SYSTEMS 

The deployment of such a large number of antennas, 
required for massive MIMO © to experience mutual or¬ 
thogonality between user channels, in confined antenna array 
dimensions results in significant increases in spatial correlation 
over conventional MIMO. In order to mitigate spatial correla- 



antenna topology for M = 256 and K = 8. 


tion, antennas can be distributed into multiple clusters, while 
assuming equal antenna array form factors. In Figures [6] and 
[7] we examine the effects of distributed antennas on linear 
precoding performance HU- 



Fig. 6: MF E[SINRj] CDF as a function of antenna cluster 
numbers, N, and antenna topology for M = 256 and K = 32. 

In Figure 0 we plot the MF E[SINR,] CDF for N = 1,2 
and 4 antenna clusters. For the i.i.d. channel, antenna cluster 
numbers have a negligible impact, other than at high SNR 
where the co-located antenna system dominates due to a larger 
number of degrees of freedom serving users in good coverage. 
When spatial correlation is introduced, larger numbers of 
antenna clusters show superior performance. This follows, 
with equal antenna form factors, as greater numbers of antenna 
clusters increase inter-element distances and thus reduce the 
effect of spatial correlation. 

Figure [7] shows the effects of varying antenna cluster 
numbers, N, on ZF E[SNR,]. Contrary to MF precoding, the 
ZF precoder has significantly better performance with larger 


































Fig. 7: ZF E[SNR,j CDF as a function of antenna cluster 
numbers, TV, and antenna topology for M = 256 and K = 32. 

antenna cluster numbers for the uncorrelated (i.i.d.) scenarios. 
This is because ^ >> K for all cases of TV considered 
and so the ZF is able to minimize all system interference 
coming from the large number of co-scheduled users in the 
system (K = 32). The gain in median i.i.d. E[SNR/ of 
approximately 20 dB, from TV = 1 to TV = 4, is therefore 
coming from the better coverage which distributed antennas 
provides, increasing the link gains, /?„,&. On the other hand, 
the MF precoders optimality is dependent on the number of 
degrees of freedom. By distributing the antennas into multiple 
clusters, we are reducing the number of effective degrees of 
freedom, e.g., a user close to the co-located antenna system 
is receiving a strong desired signal from M = 256 antennas, 
whereas a user close to the TV = 4 antenna clusters system is 
only receiving a strong desired signal from M = 64 antennas. 
This decrease in effective number of degrees of freedom 
reduces the ability of the MF precoder to mitigate inter-user 
interference, which manifests itself as off-diagonal elements 
in (G t G *)/M. Note, that in the limit of an infinite number 
of degrees of freedom, user channels become orthogonal, i.e., 
inter-user interference is eliminated, and the the performance 
of the MF precoder will approach ZF. 

V. Massive MIMO Channel Model 


accounts for NLOS propagation and H^^los) = l« x i is an 
4 x 1 vector of ones which accounts for LOS propagation. 
Kf is the Rician K-factor which controls the dominance of the 
LOS component relative to the NLOS component. 



Fig. 8 : MF E[SINRj] CDF as a function of Rician K-factor, 
Kf, and antenna topology for M = 256 and K = 8 . 



and antenna topology for M = 256 and K = 8 . 


With the deployment of small cell technology coexisting 
with massive MIMO in next generation wireless systems, there 
is a higher probability of LOS transmission lH6l . The effects 
of LOS propagation on large antenna array systems has not 
been well studied and thus we explore its impact on massive 
MIMO system performance. We model the presence of LOS 
via a Rician channel with varying K-factor. 

The j x 1 i.i.d. Rician channel vector H n ,k, between the 
/ 7 th antenna cluster and the fcth user, is given by |fI31 , IfTTl 


Hn.I — 


1 


K f + l 


-H 


ra.fc(NLOS) 


K f 


1 + Kf 


-H 


n,k( LOS)? 


(17) 


where the ^ x 1 vector FR^nlos), with CAf( 0,1) entries, 


The impact of varying amounts of LOS propagation on 
MF E[SINR,] and ZF E[SNR,:] is shown in Figures [ 8 ] and 
[9] respectively. For both antenna topologies the M = 256 
transmit antennas are configured as 16 x-pol antennas in the 
^-dimension and 8 x-pol antennas in x, y-plane. For the URA, 
this corresponds to a A /8 and A/4 adjacent inter-element spac- 
ings respectively. The small inter-element spacings produce a 
large immediate off-diagonal entry in R, given in (fl2l) . of 
nearly 0.94. 

It is observed in Figure [ 8 ] that an increased LOS presence, 
or increasing K-factor, drastically reduces the MF E[SINR.;] 
performance. For example, increasing the K-factor from Kf = 
0 to Kf = 2 decreases the median MF E[SINR*] of an i.i.d. 




























channel by approximately 22 dB. Even with a very large inter¬ 
element spatial correlation factor at 0.94, the effects of any 
LOS propagation show a much more detrimental effect over 
spatial correlation for MF E[SINR,] performance. 

As compared with MF precoding, in Figure[9]it can be seen 
that ZF is more resilient to the effects of LOS propagation. The 
presence of LOS forms another kind of correlation which the 
ZF precoder is able to minimize, thus reducing its detrimental 
effect. However, any LOS propagation conditions produce sub- 
optimal performance. 


[16] C. T. Neil, M. Shall, P. J. Smith, and P. A. Dmochowski, “Massive 
MIMO convergence analysis for next generation wireless systems,” IEEE 
Communications Magazine., to be submitted. 

[17] C. Xiao, Y. R. Zheng, and N. C. Beaulieu, “Novel sum-of-sinusoids 
simulation models for Rayleigh and Rician fading channels,” IEEE 
Transactions on Wireless Communications, vol. 5, no. 12, pp. 3667- 
3679, December 2006. 


VI. Conclusion 

In this paper, we have shown the impact of various practical 
aspects on massive MIMO system performance. We found 
that multiple clusters reduces the effects of spatial correlation, 
and in turn improve linear precoding performance. There is 
shown to be a negligible difference in system performance 
between the cylindrical and URA antenna topologies. Also, 
the ZF precoder is able to minimize the virtual interference 
introduced by LOS presence and spatial correlation, resulting 
in superior performance over MF precoding. 
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